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We report a degenerate quasicrystal in Monte Carlo simulations of hard triangular bipyramids 
each composed of two regular tetrahedra sharing a single face. The dodecagonal quasicrystal is 
similar to that 1 recently reported for hard tetrahedra [Haji-Akbari et al., Nature (London) 462, 
773 (2009)] but degenerate in the pairing of tetrahedra, and self-assembles at packing fractions 
above 54%. This notion of degeneracy differs from the degeneracy of a quasiperiodic random tiling 
arising through phason flips. Free energy calculations show that a triclinic crystal is preferred at 
high packing fractions. 



Hard disks and spheres order into hexagonal and face- 
centered cubic crystals, respectively, above a certain 
packing fraction. A more complex phase behavior is ob- 
served if the disks or spheres are rigidly bonded into 
dimers (dumbbells) [1-4]. A solid phase, disordered in 
the orientation of dimers while ordered on the monomer 
level, forms if the distance between monomers within a 
dimer is roughly the diameter of a monomer. This equi- 
librium solid phase can be alternatively understood as a 
random pairing of neighboring monomers within the na- 
tive monomer crystal. The resulting thermodynamic en- 
semble of ground states is degenerate and the structure is 
therefore called a degenerate crystal. As shown by Woj- 
ciechowski et al. [1] for hard disks, the entropy associated 
with the degeneracy exceeds the entropy from excluded 
volume effects, which by itself is sufficient to drive the 
crystallization of hard monomers. Other consequences of 
the pairing of monomers into dimers include topological 
defects [5], a restricted, glassy dislocation motion [6, 7], 
and unusual elastic properties [8]. Similar degenerate 
phases have also been observed for freely-joined chains of 
hard spheres [9, 10]. 

Although degenerate crystals can potentially assem- 
ble from dimers of hard shapes other than disks and 
spheres, few examples have been reported. One reason 
is the competition between degenerate crystals and the 
liquid crystalline phases frequently observed for particles 
with large aspect ratios. For example, elongated tetrag- 
onal parallelepipeds, which for an aspect ratio of 2:1 can 
be viewed as dimers of face-sharing cubes, form a de- 
generate parquet phase at intermediate densities before 
transforming into a smectic liquid crystal that eventually 
crystallizes [11]. Another simple dimer is the triangu- 
lar bipyramid (TBP), which consists of two face-sharing, 
regular tetrahedra (Fig. la). The TBP is the simplest 
face-transitive bipyramid and the twelfth of the 92 John- 
son solids. The lack of inversion symmetry of the TBP, 
however, makes lattice packings non-optimal [12], and 
thus it is potentially more interesting as a dimer than 
dimers of spheres and cubes. Moreover, the recent syn- 
thesis of TBP-shaped nanoparticles and colloids [13-16] 
makes the investigation of this building block of practical 




FIG. 1: Phases formed by (a) triangular bipyramids (TBPs): 
(b) TBP crystal, (c) degenerate quasicrystal, (d) regular qua- 
sicrystal approximant, (e) degenerate quasicrystal approxi- 
mant. For visualization purposes, we show member tetrahe- 
dra of most TBPs at 30% actual size and connect their cen- 
ters with bonds. In (c-e), tetrahedra and bonds are colored 
according to their orientation projected on the plane. 

relevance. 

In both of the known ordered phases of hard, reg- 
ular tetrahedra, each tetrahedron is in almost-perfect 
face-to- face contact with at least one other tetrahedron. 
The densest known packing of tetrahedra (</> = |^ ~ 
85.63%) is a parallel arrangement of two dimers (four 
tetrahedra) - that is, two TBPs - in a triclinic unit cell 
to form a dimer crystal [17, 18], which we refer to in the 
present paper as the TBP crystal (Fig. lb). At lower 
packing fractions, hard tetrahedra assemble into a do- 
decagonal quasicrystal [19], in which the tetrahedra form 
a decorated square-triangle tiling [20] . Degenerate phases 
are impossible in the TBP crystal because the contacts 
between neighboring tetrahedra in different TBPs are 
highly imperfect [21], but are possible in the quasicrystal 
due to the almost-perfect face-to-face contacts between 
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FIG. 2: (a) TBPs assemble into a dodecagonal quasicrys- 
tal in isobaric and isochoric Monte Carlo simulations, (b) 
The square-triangle tiling obtained by connecting the cen- 
ters of 12-fold rings of member tetrahedra. Intra- and inter- 
TBP bonds are depicted in black and gray respectively. (d,e) 
Diffraction patterns with centers of member tetrahedra as 
scatterers calculated (c) perpendicular to and (d) across the 
layers, (e) Intra-TBP and (f) total bond order diagrams. 

all neighboring tetrahedra. Quasicrystals are solids with 
long-range order but without periodicity [22]. Originally 
discovered in metallic alloys [23], many alloy quasicrys- 
tals are now known, and a handful of quasicrystals have 
been reported in non- metallic systems. Among them are 
quasicrystals made from spherical micelles [24], binary 
nanoparticles [25], and hard tetrahedra [19]. 

In this Letter, we investigate the phase behavior of 
hard TBPs and report a degenerate quasicrystal. The 
notion of degeneracy discussed here should not be con- 
fused with the extensively studied degeneracy associated 
with random tiling quasicrystals [26, 27] where tiles with 
unique decoration patterns mix to form random tilings. 
We instead report a new type of randomness in the level 
of decorating individual tiles, in addition to the degener- 
acy of the random tiling. 

We use isochoric and isobaric Monte Carlo (MC) sim- 
ulations to study hard TBPs, which we model as perfect 
polyhedra with sharp vertices and edges of unit length a. 
Simulations are carried out within periodic boxes with 
system sizes ranging from 432 to 8,000 particles. Each 
isochoric MC cycle comprises one update per particle on 
average, which is either a trial translation or a trial ro- 
tation with equal probabilities. An additional box trial 
move is included per isobaric cycle. For fluid phases, 
the box is resized isotropically only, while for crystals 
its shape is also allowed to fluctuate. Free energies are 
calculated using thermodynamic integration and a modi- 
fied Frenkel-Ladd method [2, 28] as described in detail in 
[21]. Further details and simulation parameters are given 
in Ref. [33]. 

The dodecagonal quasicrystal of TBPs forms sponta- 



neously from the equilibrium fluid phase at packing frac- 
tions above 54%. Fig. 2a depicts a side view of the 
quasicrystal formed in an isobaric simulation of 2,624 
TBPs at reduced pressure P* =Pa 3 /k B T = 46 and sub- 
sequently compressed to a packing fraction of 81.34%. 
TBPs arrange into layers (white lines), which stack on 
top of each other perpendicular to the 12-fold symmetry 
axis (dark arrow). We confirmed that the formation of 
the quasicrystal occurs reproducibly in systems with at 
least a few thousand particles and does not depend on 
the shape of the simulation box. 

The quasicrystal structure can be best understood by 
replacing each bipyramid by its two member tetrahedra. 
Fig. 2b depicts the centroids of tetrahedra within a few 
layers of Fig. 2a. Neighboring tetrahedra are connected 
with bonds [34] . Dodecagons that are depicted in purple 
in Fig. 2b correspond to rings of twelve member tetra- 
hedra, a structural motif characteristic of the quasicrys- 
tal [19]. These rings are further capped with pentagonal 
dipyramids (PDs), five tetrahedra sharing an edge, visi- 
ble in the figure as pentagons (green) within dodecagons. 
Additional member tetrahedra, referred to as intersti- 
tials, fill the space between the rings and are depicted in 
dark blue. Together, dodecagons and PDs form motifs 
whose centers are the vertices of square and triangle tiles. 
Their mixing gives the square-triangle tiling its overall 
twelve-fold symmetry as observed in the diffraction pat- 
tern depicted in Fig. 2c. Layering along the 12-fold axis 
can be seen in Fig. 2d. Overall, the arrangement of the 
member tetrahedra is identical to that reported in the 
hard tetrahedron system [19]. 

To elucidate how the bipyramids are arranged within 
the quasicrystal, we compare statistical distributions of 
intra-TBP bonds (bonds that connect member tetrahe- 
dra within TBPs) and the set of all bonds in the qua- 
sicrystal by projecting both sets onto the surface of a 
unit sphere. The resulting diagrams are referred to as 
intra-TBP and total bond order diagrams, respectively 
and are visualized using the Mollweide projection with 
the 12-fold axis pointing in the vertical direction. Com- 
paring these bond order diagrams (Figs. 2e,f), we observe 
no significant difference in the distribution of bond di- 
rections within the 12-fold layers. This suggests that the 
pairing of tetrahedra in the quasicrystal does not follow 
a predefined set of rules and is instead random. How- 
ever, tetrahedra tend to pair more strongly within layers 
than between neighboring layers, a fact that can be ex- 
plained by noting that face-to-face contacts are more per- 
fect within layers. Motivated by studies of hard sphere 
dimers [1], we refer to the TBP quasicrystal as a degen- 
erate quasicrystal (DQC). The randomness can be seen 
clearly in Fig. lc. It is surprising that the structural 
quality of the DQC is uncompromised despite the addi- 
tional geometrical constraints imposed on the system by 
pairing tetrahedra into TBPs. For instance, we find that 
the maximum packing fraction achieved by replacing the 
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FIG. 3: (a) Top and (b) bottom views of the regular ap- 
proximant. The unit cell has 41 triangular bipyramids. 
Particles are colored according to their environment: PD- 
PD (light green), PD-ring (dark green), intra-ring (purple), 
ring-ring (red), ring-interstitial (blue), interstitial-interstitial 
(cyan), (c) Schematics of the unit cell with connections be- 
tween neighboring rings and between rings and central in- 
terstitials shown with red and blue double- arrows, respec- 
tively. (d,f) Intra-TBP and (e,g) total bond order diagrams 
for (d,e) the regular approximant and (f,g) the degenerate 
approximant. In the legends, 'hit' stands for 'interstitial'. 

bipyramids with individual member tetrahedra and then 
compressing is statistically identical to that obtained in 
simulations of hard tetrahedra. 

Approximants are periodic phases that are structurally 
similar to the quasicrystal locally [29]. Constructing an 
approximant of the TBP quasicrystal involves not only 
choosing a periodic tiling and decorating it with tetra- 
hedra, but also pairing the tetrahedra into bipyramids. 
We choose the (3.4.3 2 .4) Archimedean tiling which, in 
the case of hard tetrahedra, gives rise to the densest ap- 
proximant [19]. There is no unique way of pairing tetra- 
hedra into TBPs even within a single unit cell of the 
approximant due to degeneracies associated with rota- 
tions of the capping PDs. In particular, it is not possible 
to avoid breaking the four-fold symmetry of the approx- 
imant unit cell in the pairing process. We constructed a 
regular approximant by retaining as much of the symme- 
try as possible. Top and bottom views of the constructed 
approximant are depicted in Figs. 3a, b while a unit cell 
is depicted in Fig. 3c where ring-ring and ring-interstitial 
connections are highlighted. We find that the regular 
approximant can be compressed to a maximum packing 
fraction of 83.39%, a bit less than the maximum packing 
fraction of 85.03% achieved for the quasicrystal approxi- 
mant constructed of individual tetrahedra [19]. The dis- 
tinctive difference between the intra-TBP (Fig. 3d) and 
the total bond order diagrams (Fig. 3e) is a result of this 
deterministic pairing (Fig. Id). 

By expanding the regular approximant, we find that 



it melts at P* < 35 and packing fractions <p < 54%. 
But before melting, the crystal slowly transforms into a 
more loosely packed structure in which tetrahedra are 
paired at random into TBPs, just as in the DQC, al- 
though their positions and orientations are unchanged 
(Fig. le). The resulting structure is therefore degenerate 
to the tetrahedron-based approximant and we refer to 
it as a degenerate approximant (DA). The angular dis- 
tribution of intra-TBP bonds around the four-fold axis 
(Fig. 3f] is more similar to that of all bonds (Fig. 3g] 
in the degenerate approximant than in the case of the 
regular approximant [Figs. 3d,e], which again suggests 
random pairing. We find that the transformation from 
regular to degenerate approximant is irreversible on the 
time scale of our simulations (~ 10 8 MC cycles). Since 
the DA can only be recompressed to a density of 82.88%, 
which is lower than the maximum density of the regular 
approximant, the DA has to be stabilized by its pairing 
disorder close to melting. 

To understand how the regular approximant trans- 
forms into the DA, we note (Fig. 2b) that the arrange- 
ment of the member tetrahedra can be alternatively 
understood as a spanning network of interpenetrating 
PDs [19]. In the hard tetrahedron system, PDs can easily 
rotate around their principal axes [21]. Such rotations 
are also essential in understanding the local rearrange- 
ments of bipyramids at densities below 60%. As shown 
in Fig. 4b, TBPs move very little at <\> = 60%. Even af- 
ter 250 million MC cycles only a small fraction of TBPs 
have moved as much as a. A much faster dynamics oc- 
curs at <j) = 57%. Particles at or near that density move 
over discrete distances that are characteristic of a PD 
network (Fig. 4a). These rearrangements change neither 
the tiling nor its decoration. Instead, they reshuffle the 
pairing pattern by a sequence of PD rotations. After 
a sufficiently large number of reshuffling moves the DA 
emerges from the regular approximant. 

Next we study the relative thermodynamic stability 
of various phases. We first compare the DQC obtained 
in simulation and its constructed approximants. As ob- 
served in Fig. 5a, both the regular and the degenerate 
approximant are slightly denser than the DQC at all 
pressures. The relation G(P 2 *) - G(Pf) oc ^dP* 
between the free energy and the equation of state then 
suggests that the approximants are thermodynamically 
preferred over the DQC at sufficiently high pressures be- 
cause their Gibbs free energies increase more slowly with 
pressure. Furthermore, the approximants melt at lower 
pressures than the quasicrystal, which indicates that they 
might even be more stable than the quasicrystal at all 
pressures. Nevertheless, the DQC remains the only or- 
dered phase that forms in our simulations. It is also the 
only structure we expect to be observed in experiments 
of hard nanocolloidal TBPs since the kinetic process of 
transforming from the DQC into the approximant is ex- 
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FIG. 5: (a) Equation of state for the TBP crystal, the degen- 
erate quasicrystal, the regular and the degenerate approx- 
imants. (b) The free energy difference between the TBP 
crystal and the approximant. (c) Equilibrium phases of hard 
TBPs. 



FIG. 4: The self part of the van Hove correlation function 
G s (r,t) measures the particle motion in the approximant. 
The separation distance r(t) is calculated between centers of 
mass of member tetrahedra. (a) Large rearrangements occur 
at (p = 57%. (b) There is little motion present at cj> = 60%. 
The observed dynamics is similar to that observed in the hard 
tetrahedron system [21]. 

tremely slow. Considering the local structural similarity 
of the DQC and the fluid in terms of the PD network, 
the formation of the less stable DQC and not the approxi- 
mant in simulation may be another example of Ostwald's 
rule [30]. 

Next, we compare the approximant with the TBP 
crystal by calculating the free energy difference between 
them. As shown in Fig. 5b, the approximant has a lower 
free energy than the TBP crystal for packing fractions be- 
low 79%. A phase transition occurs at P* = 356 ±50, cor- 
responding to coexistence packing fractions of (/> c ,app = 
(79.1±0.8)%and^ c ,TBP = (80.7±0.7)%. The thermody- 
namic stability of the approximant at lower densities can 
be attributed to the additional configurational entropy 
associated with collective motions of particles. Such mo- 
tions are not present in the TBP crystal. Their role in 
stabilizing the quasicrystal approximant has been shown 
for the structurally and dynamically similar system of 
hard tetrahedra [21]. The phase diagram of the hard 
TBP system is depicted in Fig. 5c. 

Remarkably, hard TBPs not only prefer a complex qua- 
sicrystal over the simpler TBP crystal at intermediate 
packing fractions, but also form it on timescales com- 
parable to that previously observed in the hard tetrahe- 
dron system. This is surprising because, in comparison to 
tetrahedra, the motion of the highly anisotropic bipyra- 



mids is considerably more constrained. Nevertheless, the 
degeneracy of the quasicrystal helps it form easily in sim- 
ulation. Random pairing allows TBPs to join existing 
seeds of the DQC without forming configurations that 
are kinetically trapped due to incorrect pairing. Parti- 
cle rearrangements needed for the formation and growth 
of the seed are also feasible due to the local similarity of 
the fluid and the quasicrystal [19, 31]. Finally, the degen- 
eracy and the existence of ring-ring and ring-interstitial 
"cross-links" adds rigidity to the TBP structures. This 
means that the TBP system might be superior over the 
tetrahedron system in terms of its mechanical proper- 
ties, just as for crystals of hard sphere dimers compared 
to crystals of their monomers [8]. 

In conclusion we have shown that hard triangular 
bipyramids form a degenerate dodecagonal quasicrystal. 
Our finding is only the second quasicrystal formed with 
hard particles, the first reported degenerate quasicrys- 
tal, and one of only a few quasicrystals formed in non- 
atomistic systems. Our results suggest that degenerate 
phases are not restricted to simple close-packed crystals 
and might be common in dimer systems. 
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SUPPLEMENTARY INFORMATION 
Compression Algorithm 

To overcome the sluggishness of conventional isobaric 
Monte Carlo (MC) simulations in compressing hard par- 
ticle systems with large numbers of particles, we utilize a 
modified compression algorithm described in detail in [19] 
for obtaining the densest packings of the quasicrystal and 
the approximants. This method relies on allowing a small 
number of minor overlaps in compression moves, which 
are subsequently removed to obtain dense packings. All 
expansion/compression moves are accepted even if they 
result in new overlaps and a different criterion is used to 
ensure that the number and amount of overlaps remains 
small. We monitor p, the fraction of trial translations 
accepted since the previous volume move, and compare 
it to a target acceptance probability p t . If p < pt, the 
volume of the box is increased by a random factor / uni- 
formly chosen from the interval [1, 1 + 0.002Ax] while for 
p > p t a compression is attempted by a random factor 
uniformly distributed in [1 — 0.002Ax,l]. Between vol- 
ume moves, the system evolves through conventional trial 
translations and trial rotations and the moves that do not 
generate new overlaps are accepted. Ax, the maximum 
step size for a trial translation is the control parame- 
ter in this method and is inversely related to pressure in 
the conventional isobaric MC algorithm. In a compres- 
sion run, Ax is exponentially decreased until the densest 
packing is obtained. A few hundred cycles of conven- 
tional NVT MC are carried out at the end to remove any 
overlaps generated by this method. 

Free energy calculations 

We use a modified version of the Frenkel-Ladd ther- 
modynamic integration scheme [2, 28] to calculate the 
Helmholtz free energy differences between various crys- 
tals in the system. The Gibbs free energy is then deter- 
mined from the Helmholtz free energy together with the 
equation of state. The equation of state is also used to 
extrapolate the Gibbs free energies to pressures where no 
Frenkel-Ladd calculation can be performed by integrat- 
ing dG = VdP. The values of 7 max and c chosen in this 
study are included in Table SI. Further technical details 
of the free energy calculation scheme can be found in [21], 
where the same method was applied to a system of hard 
tetrahedra. 

Van Hove correlation function and structure factor 

The van Hove correlation function and the structure 
factor are calculated for centroids of member tetrahedra 
in the hard TBP quasicrystal approximant. To deter- 



mine the structure factor, the centroids are convoluted 
with a Gaussian and projected along the observation di- 
rection. The resulting pattern is sheared into a square, 
discretized, fast Fourier transformed, and then sheared 
back. 



Correlated Motions and Additional Entropy 

The entropy of a hard particle system is the volume 
of the 6Af-dimensional configuration space accessible to 
it. In dense arrangements where particles rattle in their 
respective cages, this accessible volume can be approx- 
imated by the product of free volumes of individual 
particles with the free volume being the volume of the 
6-dimensional configurational space of a given particle 
accessible to it while all other particles are kept fixed 
(mean-field approximation [21, 32]). Through collec- 
tive motions, particles can access regions of the phase 
space not accessible through rattling of individual par- 
ticles. Therefore collective motions contribute to some 
additional entropy not accounted for by the mean-field 
approximation and further stabilize the corresponding 
system. The stabilizing effect of correlated motions in 
the quasicrystal approximate was rigorously confirmed 
for the hard tetrahedron system by comparing exact free 
energies calculated from the Frenkel-Ladd method with 
the free energies estimated from calculated free volumes 
and the mean-field approximation [21]. A comparison of 
Fig. 4 in the present work with Fig. 9 in [21] demon- 
strates the remarkable similarity of the dynamics in the 
TBP and tetrahedron quasicrystal approximants, in par- 
ticular with regards to collective motions of the type de- 
scribed above. We therefore expect these collective mo- 
tions to have similar stabilizing effects in both systems. 

Dense packings 

The densest unit cells of the TBP crystal and the reg- 
ular approximant are given in Tables S2 and S3, respec- 
tively. In the tables, the orientation of each TBP is rep- 
resented by a unit quaternion q = (q t: q X: q y , q z ) with 
rotation matrix 

q 2 + q 2 x -q 2 y -q 2 z ^(QxQy-qtqz) 2(q x qz + qtq y ) \ 
2(q x q y + qtq z ) q 2 - q 2 x + q 2 y - q 2 z z(q y q z -qtq x ) 
2(q x q z - qtq y ) z(q y q z + qtq x ) q 2 -q 2 x -q 2 y + q 2 z ) 

The vertices of the TBP with q = (1,0,0,0) are given 
by: 

vi,2 = (o,0,±4V3/3),v 3) 4=(-V6/3,±V^,o), 
v 5 = (2v / 6/3,0,o). 
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TABLE SI: Simulation details for the calculation of the equation of state (EOS), quasicrystal (QC) assembly, and thermody- 
namic integration. Several independent runs were performed for each state point to assure accurate statistics. The smallest 
system used for quasicrystal formation has 1,458 particles. Quasicrystal formation is robust and is routinely observed at densi- 
ties above 54% and system sizes larger than 1,400 particles. We do not observe a quantitative difference in the structure of the 
quasicrystal for such large system sizes. The small, 432-particle system is only used for the mathematically constructed TBP 
crystal (with 216 two-particle unit cells), which is used for estimating the equation of state and free energy calculations only. 

Parameters 
0.01 < P* < 60 
40 < P* < 60 
0.5 < < 0.6 
32 < P* < 10,000 
36 < P* < 10,000 
35 < P* < 10,000 
0.57 < cf) < 0.80 
< 7 < 4 x 10 6 ,c= 1/2 

0.70 < < 0.80 
< 7 < 4 x 10 6 ,c= 1/2 
0.60 < 6 < 0.80 



Phase Objective 


Ensemble System size 


MC sweeps 


Fluid EOS calculation 


isobaric 


2,624 


10 7 


Fluid QC assembly 


isobaric 1,458 - 2,624 


10 8 


Fluid QC assembly 


isochoric 


8,000 


10 8 


TBP crystal EOS calculation 


isobaric 


432 


10 7 


Quasicrystal EOS calculation 


isobaric 


2,624 


10 7 


ApproximantEOS calculation 


isobaric 


656 


10 8 


ApproximantDynamics 


isochoric 


656 


3 x 10 8 


TBP crystal Therm. Integration isochoric 


432 


2 x 10 5 per 7 


Approximant Therm. Integration isochoric 


656 


2 x 10 5 per 7 



TABLE S2: A unit cell of the TBP crystal with = 85.6347%. The lattice vectors are given by 61 = (8/5, 12/5, 4/5), 62 = 
(3/8, 97/80, 191/80), 63 = (12/5, 3/20, 3/2). The positions and quaternions of the two particles in the unit cell are given. 



i 


Xi 


Vi 




Qt,i 


Qx,i q y ,i 


Qz,i 


1 











0.880476- 


-0.3647050.279848- 


-0.115917 


2 


-26/15- 


-79/120- 


1/120 


-0.279848 


0.1159170.880476 


-0.364705 



8 



TABLE S3: A unit cell of the regular approximant with = 83.39%. Lattice vectors are given by bi = (9.96864731, 0, 0), b 2 = 
(0, 10.12720286, 0), b 3 = (0, 0, 2.597424124). There are 41 particles in the unit cell; their positions and orientations are given. 
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—4 


"1 "1 A A A A A r\ 

118404870- 


-0 


A A A A O /° 1 -1 A 

892236112 





A A f A A 7" /"» A 
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A i"\ C\ A O r" A A /I 
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-0 


A A A ~1 A r* r* 
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<A A A A O A r" A 
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OAOAAf" A P* P~ 

393025455 


3 


-3 


924682848- 


-3 


466017697- 


-0 


t* a a 1 -1 r* 

350811567 




-0 


O-l AA-1 1*7 A A /"» 

318217496 





a 6^% a a r\ a a a 

434494644 





f7A rn a r* a 1*7 p* 

725846175 





A A f*7 A A1 111 

42790HH 


4 


— 


A f C\A OO /I /~' O 

462122468- 


-4 


240376441- 


-0 


782025752 
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387765410- 


-0 


303556385 





174836762 


5 


— 


O A A A C^t A a r* 

327948406- 


-1 
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-1 
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f*7 r* A A A f*7 f*7 f*7 O 
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—2 
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-0 
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7 


—2 


069478657- 


-1 


471298448- 


-0 


432573444 







-1 T^/^r aoai 

177650291 





552837588 





A A f * 1 A /I A O A 

096194039 





O AO /I O A 1 A f 

808429145 
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A A Z'"* A /I A A y| 

826248624- 


-0 


C\ 1*7 1*7 A A 7" f" A A 
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-0 
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-4 
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-0 
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O 1 A 1 /I O A /I /I 
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A A O A Af7 A f O 

142097453 


11 


—4 


69664H86 
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a^7/° aa 7" r* a 
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/^*AAA /I AA^70 

602240973 







O r* A A A A A A A 

360108299 





A f" A K A 01*71*7 A 

252500774 





A/^*A /I AO^TAA 

269483792 





Ar/^^ArAO yl 

856705234 


12 


—3 


747047319 
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Ar 7" A A 1 A /I A 

055281040- 


-1 
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^TAAAAAf" A^T 

780999587- 


-0 
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449388327 





430045356- 


-0 


Ar/^*1 OI l*7 A A 

056131712 


13 


-0 


058997540 


4 
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025246097- 


-0 
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-0 
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14 
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019924758 
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-0 


681742157 




-0 
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AAA/"' O 1 /I O 
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15 
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691606294 





561167807- 


-0 


315127681- 


-0 


327820027 
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431835202 
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22 


4 


166061318- 


-1 


^rr" 7" a a /i 

755893743 





r* A A /I A A r" A A 

614482584 







A f ^ O O A O A 

857339372 





268730116- 


-0 


A /I OP AA A C"i A 

248620481- 


-0 


O A A r 7T A A 

361857958 
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FIG. SI: Tiling obtained for a slab of the degenerate quasicrystal formed in an NVT simulation of 8,000 particles at <j> = 54%. 
Like Fig. 2b (main text) the centers of member tetrahedra are depicted and they are colored according to their environment: 
twelvefold ring (purple), capping PDs (green) and interstitials (blue). Compared to the quasicrystal shown in Fig. 2b (main 
text), the assembled quasicrystal has more defects and is yet to fully crystallize in the layers above and below the ones shown in 
the figure. This explains the abundance of rhombs and zippers that have shown to be important in the process of crystallization. 



